variance <- function(data, na.rm = TRUE) {
	data <- as.data.frame(data)
	if (na.rm==TRUE) n <- nrow(na.omit(data)) else n <- nrow(data)
	return(sd(data,na.rm=na.rm)^2*(n-1)/n)
}

